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Non-equilibrium transport through a quantum dot with one spin-split single-particle level is studied in the 
cotunneling regime at low temperatures. The Coulomb diamond can be subdivided into parts differing in at 
least one of two respects: what kind of tunneling processes (i) determine the single-particle occupations and (ii) 
mainly contribute to the current. No finite systematic perturbation expansion of the occupations and the current 
can be found that is valid within the entire Coulomb diamond. We therefore construct a non-systematic solution, 
which is physically correct and perturbative in the whole cotunneling regime, while smoothly crossing-over 
between the different regions. With this solution the impact of an intrinsic spin-flip relaxation on the transport 
is investigated. We focus on peaks in the differential conductance that mark the onset of cotunneling-mediated 
sequential transport. It is shown that these peaks are maximally pronounced at a relaxation roughly as fast as 
sequential tunneling. The approach as well as the presented results can be generalized to quantum dots with few 
levels. 
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I. INTRODUCTION 

Coherence effects, whose signatures can be seen in (spin-) 
electronic transport through low dimensional nanoscopic 
structures like quantum dots''^'\ provide insight into funda- 
mental aspects of quantum mechanics and have important ap- 
plications in vital fields of research such as spintronics, quan- 
tum computing, and data storage."*''' For the occurrence of 
these effects spin-flip relaxation is widely considered as a lim- 
iting factor and therefore usually sought to be as small as 
possible. In experiments spin-flip relaxation times Ti rang- 
ing from /is^'^ to ms'^ '' have been observed displaying depen- 
dences both on the sort of quantum dot and on the parame- 
ters of the experimental setup like temperature and magnetic 
^gjjj2,3,io,ii Accordingly, the relaxation rates can be experi- 
mentally adjusted in a wide range either by means of tuning 
of external parameters or by suitably tailoring the quantum 
dot itself. Recently, spin-flip times of even several hundred 
milliseconds were measured in n-doped (In, Ga)As/GaAs 
quantum dots charged with spin-polarized electrons at low 
magnetic field and temperature.'^ For transport through few- 
electron quantum dots in the presence of intrinsic spin relax- 
ation, which is discussed in this paper, microscopic mecha- 
nisms like phonon-induced spin decay due to spin-orbit or hy- 
perfine interaction have been investigated theoretically (see, 
e.g., Refs. 2,10,11,13,14,15). 

Though relaxation mostly acts destructively on coherent 
electron dynamics, it can also, however, considerably pro- 
nounce their effect, as Weymann et al}^ show for the case 
of Coulomb-blocked transport through a single-level quantum 
dot (SLQD) coupled to ferromagnetic leads with antiparallel 
magnetization. The zero-bias anomaly of the differential con- 
ductance and the conductance step at the onset of inelastic 
cotunneling are increased by a slow spin relaxation for a spin- 
degenerate and a spin-split dot level, respectively. 

We will show in this paper that a similar, strongly pro- 
nounced effect should be observable even when the leads 
are non-magnetic. In the considered case a small relaxation 



that is roughly as large as the tunnel coupling maximizes 
peaks of the differential conductance, which mark the onset 
of cotunneling-mediated sequential transport. This effect is 
associated with sequential tunneling out of an excited single- 
particle state. Within a SLQD model, the dot's level has to 
be spin-split. For few-electron GaAs/Alo.aGao.yAs quantum 
dots with non-degenerate orbital levels these signatures of the 
single-particle spectrum have been intensely studied both ex- 
perimentally and theoretically by SCHLESER et al. whereby 
new insight was provided into the interplay between sequen- 
tial and cotunneling in the Coulomb blockade regime. In the 
present paper we investigate this interplay in further detail. 
Since said signatures appear close to resonance with single- 
particle transitions, we have to base our calculations on the 
non-equilibrium Keldysh formalism rather than second-order 
perturbation theory. '^"'^"'^ It is shown that in order to obtain 
physically correct, perturbative results for the entire cotunnel- 
ing regime, one has to construct non-systematic rate equations 
similar to those proposed in Ref. 20. In the latter equations 
we identify terms that cannot belong to the second-order per- 
turbation expansion. By omitting these terms, one ensures 
that for the considered system the occupation probabilities are 
well-defined everywhere in the Coulomb blockade regime. As 
in Ref. 16 we treat the effect of relaxation phenomenologi- 
cally, describing the intrinsic spin-flip processes by an effec- 
tive rate 6. Thus no particular mechanism has been specified. 

The paper is structured as follows. In Sec. II we introduce 
the model system and explain restrictions on the system pa- 
rameters. The diagrammatic transport theory and the deriva- 
tion of transport equations are sketched out in Sec. III. Results 
are presented and discussed in Sec. IV and followed by a sum- 
mary in Sec. V. 



II. MODEL 

We consider a model system consisting of a SLQD, which 
is coupled to two metallic leads (L and R) by identical tun- 
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neling barriers, so that a dc bias voltage Hias, symmetri- 
cally applied between both reservoirs, causes a tunneling 
current through the dot. An additional capacitatively cou- 
pled gate electrode allows to adjust the electrostatic poten- 
tial in the dot by applying a gate voltage. Such a sys- 
tem can be represented by the Anderson-type Hamiltonian 
H = Hd + Hl + Hn + Ht with the quantum dot part 
He), the Hamiltonians H]^ and of the left and right lead, 
respectively, and the tunneling operator Ht, describing the 
coupling between the dot and the leads. We assume that 
the spin degeneracy of the two single-electron dot states is 
lifted (e.g., by a Zeeman-field), leading to an energy dif- 
ference of A. Then the dot Hamiltonian can be written as 
Hd — eagOg + (e + ^) Oefle + U a\a^a\ag. Here the index 
g (e) denotes the spin of the single-electron ground state |g) 
(excited state |e)) and a\ (0^) with a G {g,e} creates (anni- 
hilates) an electron with spin a and energy e^r — e + Sa-.eA 
when acting on the empty dot state |0) (6 is the Kronecker 
delta). U is the Coulomb-energy of the doubly occupied 
state |d). The leads play the role of macroscopic reservoirs 
and are described as free electron gases with Hamiltonian 

= J2k.cr '^k,r c\ .^Ck ,a ,t , wherc r G {L,R} refers to the 
lead; k is the wave vector of an electron in reservoir r, a is 
its spin and ek,r its energy. The c^f^ „ ^ {ck.a,r) are the cor- 
responding creation (annihilation) operators. Due to the ap- 
plied bias voltage, /i^ = (— l)'^'-''eVbia,s/2 — with e > being 
the elementary charge — gives the electrochemical potential of 
reservoir r. The coupling between the leads and the dot is de- 
scribed by Ht — J2k a ri.^ ^'a^k,a,r + h.c}j, whcrc the first 
(second) term on the rhs describes tunneling into (out of) the 
dot with the complex tunneling parameter 7 (7*), which is as- 
sumed to be independent of wave vector and spin of a tunnel- 
ing electron as well as the reservoir out of which (into which) 
it tunnels. With the constant density of states p of the reser- 
voirs the coupling can be characterized by the positive scalar 
parameter F :— p alone. ^' The stationary tunneling cur- 
rent / is the expectation value of the current operator / :— /l, 
where = -i{e/h) J2k,a{^ ^l^k,a,r - h.c.y^. 

We demand that the reservoirs stay in equilibrium even 
when coupled to the SLQD. For a perturbative calculation of 
the occupation probabilities and the current up to second or- 
der in the tunnel coupling, T has to be very small compared 
to the dot energies e and U. Coulomb blockade of sequen- 
tial transport is possible, if the thermal energy is very much 
smaller than U, i.e., f3~^ = ksT <C U with T being the tem- 
perature and /cb Boltzmann's constant. We restrict our study 
to parameter sets with PV <C 1, which is a necessary con- 
dition for physical behavior of the second-order perturbation 
expansion in T once the electro-chemical potential of a reser- 
voir is close to resonance with the energy of a single-charge 
excitation. Furthermore, to be able to see transport signa- 
tures of the excited state within the Coulomb blockade regime, 
the Zeeman-splitting A must not be very much smaller than U 
but roughly of the same order of magnitude (though not larger 
than U/2). This implies /3A ^ 1. For the particular pa- 
rameter set (r = 4.5 X IQ-^kBT, A = ibk^T, U = 22bk^T) 
we use throughout the following discussions, this requirement 



may be difficult to meet for quantum dots made of GaAs or Si 
and magnetic fields available in laboratories. On the other 
hand, the presented perturbative framework can be applied to 
systems with F that is up to 10 times larger and with [/, A that 
are 10 times smaller, while yielding the same qualitative re- 
sults. In practice, since it is purely of mathematical origin, the 
criterion (3T <C 1 is not experimentally relevant and imposes 
no restriction on the physics underlying the transports effects, 
we present here. Therefore, the results of our approach can 
also be applied to experiments on GaAs or Si quantum dots 
in which the split-exceeded level can be resolved in transport 
spectroscopy.^"* InAs nanowire quantum dots, however, have 
an effective g Factor between 8 and 9^^ and an effective mass 
of m* = 0.02me {ra^ is the mass of a free electron).-* In 
experiments with these dots an adequately large spin-splitting 
should be feasible for Zeeman fields in the range of 1 to 10 T, 
even for the parameter set we use here. 

We would also like to emphasize in this context, that the 
presented approach is not restricted to quantum dots with one 
spin-split single-particle level but can in the same way be 
employed for dots with two non-degenerate spinless orbitals. 
For very similar systems (few-electron GaAs/AlGaAs quan- 
tum dots) the discussed conductance peaks were seen in low 
magnetic fields.'^ 

III. MASTER EQUATIONS 

For our calculations we use the real-time diagrammatic 
technique developed by Schoeller et al.^^ It is based on the 
Keldysh formalism and allows to represent a dynamical, non- 
equilibrium property of the model system by a formaly exact, 
infinite perturbation expansion with small parameter F. From 
such an expression one can obtain a systematically expanded 
quantity up to a finite order in the coupling. In 
order to construct perturbative solutions for the occupation 
probabilities and the tunneling current, we first assume that 
intrinsic relaxation is absent. To compute the time-dependent 
statistical expectation value Tr(y5(t)/'") of the current opera- 
tor, we have to calculate the density matrix p{t), which con- 
tains the complete system dynamics. Since the reservoirs are 
assumed to stay in equilibrium at all times, the density ma- 
trix's reservoir degrees of freedom can be integrated out using 
Wick's theorem. This yields the reduced density matrix poit), 
which depends only on the dot degrees of freedom. Via an 
adiabatic switching between times and t the initial state of 
the isolated dot, represented by = poito), is connected to 
the reduced density matrix of the coupled system pu {t) ■ This 
relation is expressed by equation 

PD{t) = Ilit,to)pl, (1) 

where Il{t, t') is a time evolution operator describing propa- 
gation of the reduced density matrix between t' and t. The 
propagator t') can be represented as an infinite sum of 
diagrams on the Keldysh contour, each of which is decom- 
posable into parts 11° corresponding to propagation that is not 
influenced by the reservoirs and irreducible self-energy parts 
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that describe coherent dynamics governed by the tunnel cou- 
pling and allow the dot to change its state-^. With the operator 
S, which consists of all irreducible diagrams, a Dyson equa- 
tion for n can be set up leading to the kinetic equation 

ft r*^ . . (2) 

+ / dt2 dtiU'>{t,t2)^{t2,ti)pDih), 
Jto Jto 

when plugged into (1). In the limit of to — > — oo and vanishing 
adiabatic switching, the derivative of Eq. (2) with respect to t 
becomes a self-consistent conditional equation for the station- 
ary reduced density matrix p'^, provided that the SLQD will 
eventually forget its initial state p^j due to the interaction with 
the macroscopic reservoirs^"*. Since we assume diagonality of 
the initial density matrix p'^, which implicates diagonality of 
p%^^, it is convenient to replace the latter by the vector P of 
the stationary probabilities = {(j>\pf)\(l)) for the dot to be 
in state 10) with (j> e {0,g,e,d}. We then replace the tensor 
operator E with the matrix W, whose elements 

W^'^^:- r dt'E'^^CtM (3) 

are interpreted as stationary rates of quantum dot transitions 
from state | (p) to state Since the total probability has to be 
conserved, the 4x4 Matrix W has a rank of three. Therefore 
the resulting self-consistent rate equation 

WV = Q (4) 

has non-trivial solutions and, together with the normalization 
condition = 1, uniquely determines P as a vector of 

probabilities.^" ''- A similar equation for the current out of 
reservoir r into the dot can be formed, if we introduce an oper- 
ator f7 , which consists of all irreducible diagrams of E, each 
having its last internal vertex replaced by an external vertex 
stemming from T^{t — 0). With a matrix W'^ , defined in 
analogy to (3), we get 

Each of the Eqs. (4) and (5) yields an infinite system of cou- 
pled equations, if we express W, W^, P, and / as expansions 
in r and sort all terms by order. The ?i*-order occupation 
vector and current are given by 

n 

= ^ VF(''p("~') 

„ (6) 
/(«) ^_e^^(;yL(Op(n-0)^. 

1=1 

With the terms T^("' and VF*"'"^ we identify those parts of 
W and W'' , respectively, that are represented by irreducible 
diagrams with exactly n tunneling lines. Each of these lines 
connects two vertices on the Keldysh contour and represents 



the wick contraction of the corresponding reservoir operators. 
The ascending orders of P and / are then calculated iteratively, 
starting with P*^"' and where the P^"-* are normalized ac- 
cording to ^0 P^"' — (5„ 0- The first- and second-order equa- 
tions describe transport caused by sequential tunneling and 
cotunneling processes, respectively. The curve of the sequen- 
tial current against the bias voltage resembles a staircase with 
thermally broadened steps formed at bias values appropriate 
for single-charge excitations. Cotunneling further broadens 
these steps-^^ and dominates the transport behavior within the 
Coulomb blockade regime (or cotunneling regime), where the 
gate voltage is tuned to charge the SLQD with one electron, 
while the bias is too small to doubly occupy or to empty the 

jJqj 18,36,37,38 




FIG. 1: (Color online), (a) Schematic picture of the diamond-shaped 
cotunneling regime showing its subdivision into areas with differ- 
ent possible tunneling processes (hatched areas) as well as into core 
(CR), shell (SR), and intermediate region (IR), in which the occu- 
pations Pg and Pc are determined either by cotunneling, sequential 
tunneling, or both (colored (shaded) areas). E stands for elastic and 
I for inelastic cotunneling, S is for sequential tunneling, (b) Cal- 
culated charging diagram (based on Eqs. (9) and (11)) of the co- 
tunneling regime with parameters A = A^k^T, U = 22^ks,T, 
r = 4.5 X 10"^fcBT and 6 = V/2. High values of the differen- 
tial conductance G = dl /dVb\xs (in units of Go = /SFe^/ft) outside 
and at the border of the diamond are clipped by the color scale (val- 
ues near the border and in the exterior: see text). The short horizontal 
line corresponds to the range of bias values in Fig. 2. Both (a) and (b) 
can be extended to regions with opposite sign of eVbias by reflection 
with respect to the eVLias = axis. 

Depending on its strength F, the coupling of a microscopic 
system like a quantum dot to macroscopic reservoirs will 
modify the dot's behavior slightly (F ^ e. U) or drastically 
(F ^ e, U). Even in the first case^**^ and for incoherent sequen- 
tial transport'*^, the coupling may be too strong to assume that 
the dot propagates as if isolated between two tunneling events. 
In such cases and all the more when transport of highly cor- 
related electrons is considered'*'"'*^, it may still be possible to 
represent the dot as isolated but with renormalized instead of 
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the bare system parameters. Strictly speaking, this is also true 
for the system studied here and correctly accounted for by the 
diagrammatic technique.'''' However, since the difference be- 
tween renormalized and bare energies scales with F, for the 
regime we investigate (F ^ ^ e, A, U) it is so small, 
that the effect of the renormalization on the discussed trans- 
port phenomena is virtually unobservable. This can be seen, 
for example, in Fig. 1(b), where a light vertical shade within 
the Coulomb blockade regime indicates the onset of inelastic 
cotunneling as soon as |eVbias| equals the renormalized excita- 
tion energy (see below), while on the eVbias axis the position of 
the bare spin-splitting A is marked. Obviously, both positions 
do not considerably deviate from each other Hence, through- 
out the following discussions, we do not distinguish between 
the bare and renormalized quantities, although all statements 
are strictly valid only for the latter ones. 

Fig. 1(a) schematically shows the (sub-)structure of the 
cotunneling regime plotted against eVbias and the dot poten- 
tial energy. For convenience the latter is given relative to 
the energy that is needed to charge the dot with one electron 
in the ground state: e<i>^ := e^]j — e. On the one hand 
the diamond-shaped cotunneling regime (Coulomb diamond) 
breaks down into three regions differing with respect to the 
kind of tunneling processes that predominantly determine the 
occupation of the single-particle states. As we discuss below 
Pg and Pe are given by sequential tunneling in the red colored 
(medium gray) shell region (SR), by cotunneling in the green 
(dark gray) core region (CR) and by a mixture of both in the 
yellow (light gray) intermediate region (IR). 

On the other hand one can distinguish three sub-regimes 
of the Coulomb diamond with different current-driving tun- 
neling processes. The quantum dot is in state |g) (Pg « 1) 
for I eVbias I < A (dotted area) and the small finite current 
is maintained just by energy-conserving elastic cotunneling 
(E) through virtual states |0) and |d) (elastic regime). Once 
I eVbias I exceeds A, inelastic cotunneling processes can excite 
the dot into state |e) while transferring energy from the reser- 
voirs into the SLQC.'^ **^ Since each of these processes ef- 
fectively carries one electron through the dot, they cause an 
additional electron flow (h-I). When passing from the core to 
the shell part of this inelastic cotunneling regime — the corre- 
sponding areas are hatched with crossed lines and dashed hor- 
izontal lines, respectively — cotunneling-mediated sequential 
tunneling out of the excited state sets on and further increases 
the current (h-S).'^ 

Before we can study the transport in vicinity of the excited 
state resonances, it is necessary to modify the rate equations 
(6), as they prove to be unsuitable to describe the occupations 
and current in the intermediate region. As Weymann et al. 
show in Ref. 20 it is due to the breakup of the cotunneling 
regime into core, shell, and intermediate region that no sys- 
tematic second-order expansion of P or / exists, which is valid 
within the entire regime. This can be explained as follows. 

In the shell region only those sequential transitions are ener- 
getically forbidden — the corresponding rates being exponen- 
tially small — that carry the dot out of the ground state |g). 
Hence, after a finite time of propagation, the dot inevitably 
gets trapped in |g) and thereby forgets its initial state. So the 



stationary occupations are essentially determined by sequen- 
tial tunneling alone. Even with all rates W^V set to zero, the 

matrix VF'^' still has a rank of three and by Eq. (4) all pj,"'' 
are fixed except for normalization. 

In the core region sequential transitions out of both single- 
particle states are forbidden. Then classically the dot can get 
trapped either in the ground or in the excited state, so that 
the single-particle occupations depend on the initial dot state. 
Consequently, they are no longer determined by the lowest- 
but by the second-order rate equation, i.e., they are essentially 
given by cotunneling. This becomes manifest in the structure 
of W^^\ which has a rank of two, when all rates w'^S' „ „ are 

(p <— g,e 

set to zero. 

Between shell and core lies the intermediate region, 
where the system continuously changes between classical 
and cotunneling-dominated occupation, respectively. But as 
well as no matrix W^^^ can be constructed that continuously 
changes its rank, no single rate equation exists that both de- 
termines the systematic second-order expansion of P in terms 
of W''^\ Ty*-^^ and is valid within all three regions simulta- 
neously. 

Alternatively, we seek second-order approximations of P 
and /, which are valid in the cotunneling regime and perturba- 
tive in the sense that they deviate from the systematic expan- 
sions at most by terms quadratic and cubic in F, respectively. 
With the normalized solution P' of Eq. (4), in which W is 
replaced by the sum of the lowest two orders VF'^^^-*, Wey- 
mann et al. present an example for an approximation that 
is perturbative even for arbitrary values of eVbias and e^fj.^" 
Unfortunately, for our system P' isn't well-defined within the 
entire core region, where the component P^. becomes nega- 
tive when I eVbias I < A (green (dark gray), dotted area in Fig. 
1(a)). To resolve this problem we take into account that on the 
rhs of the second-order equation'*"' 

p(A(i)p(i) = _VF^^^P("' (7) 

the first-order probabilities Pq*'], which are exponentially 
small within the cotunneling regime, are multiplied with 
the rates VF^,^q j. Hence, these rates drop out of Eq. 
(7) and its rhs reduces to a vector V with components 

V^' = - J2^i^<i>,g + ^0,e)W^0'<_0- As a consequence the 

(2) 

rates W^,1_q ^ don't contribute to systematic expansion orders 

given solely in terms of VF*-^' and W^'^\ and all contribu- 
tions to P', they are contained in, are unsystematic and should 
be omitted. With regard to the approximation of the current 

h(2) 

the same is true for terms containing the rates W^, J_q ^. By 
dropping the unsystematic terms we arrive at 

= (VT^^) + Ty^^^)P, (8) 
/=-e5]((W^L« + l^"^'^)P)^, (9) 

where W^^l^ = (S^^, + <5^,e)VF(l , . 
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Finally, without specifying a particular spin-flip mecha- 
nism, we include relaxation via an effective Hamiltonian 



0.2 0.4 0.6 0.8 1 



i-el 



h.c), 



(10) 



which describes the coupling of the dot electrons to a bath of 

free particles with temperature T and i/bath — J2q^q^q^q'^^ ■ 
This coupling is characterized by the single complex parame- 
ter T, giving the amplitude for a spin-flip process from |e) to 
|g). We assume, that the relaxation processes are completely 
incoherent to the electron tunneling and include only the first 
order of the perturbation expansion with respect to H^d- Then, 
in the diagrammatic representation, the self-energy operator 
up to second order becomes the sum of all irreducible dia- 
grams that have either up to two tunneling lines (S'^'-^') or 

exactly one relaxation line (S^gj^), which represents a wick 
contraction of bath operators. The latter operator gives rise 
to an additional matrix term in the master equation, whose 
matrix elements are defined in analogy to Eq. (3). Hence, we 
get the rate equation 



0= 



0)P 



(11) 



for a relaxation-dependent approximation P. Evaluation of the 
relaxation diagrams then yields the rates 



©g^ — e — 0e< — e 



Ge^g — Qg^g 



27r|T| 

h 

2ti\t\ 



d<^q{bqK)bPb{eq)S{eq - A) 

deq{blbq)hPb{£q)S{eq - A), 
(12) 



where Ph{^q) gives the density of states in the bath at energy 
Eg and (•)b denotes the expectation value with respect to the 
bath degrees of freedom. Alternatively, these spin-flip rates 
can be calculated using standard time-dependent perturbation 
theory and Fermi's Golden rule (see, e.g., Ref. 46). The first 
equalities in the Eqs. (12) express the conservation of the total 
probability, which in the diagrammatic approach is fulfilled by 
construction. Since we assume that the spin-splitting is large 
compared to the temperature (A/3 ^1), the relaxation rates 
are approximately given by Q^'^^ = S^A^^'.g ~ S^,^^)9/h 
with 6 = 2tt\t\^ Pb(A) both for a fermionic and a bosonic 
bath (as long as |/ibath| ^ A)."*^ P as well as / are well- 
defined and perturbative within the cotunneling regime and 
seamlessly link in the intermediate region the systematic ex- 
pansions that are only valid either in the core or shell. We also 
note that, because iJrei and Hr commute, the relaxation does 
not contribute (directly) to the current /, that is to say, Eq. (9) 
remains valid without modification. 



IV. RESULTS 

In this section we argue that the rich internal structure of 
the Coulomb diamond with its different overlapping regions 
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FIG. 2: (Color online). Current I (in units of Jo = eV/h) and 
differential conductance G versus bias eVbias and relaxation 9 with 
e "3?^ = ([/ + A) /2 and parameters A, U, F as in Fig. 1. It can be 
seen how the height of the cotunneling-mediated current step (a) and 
of the corresponding conductance peak (b) as well as their positions 
depend on the relaxation rate. In the range of < 6 < F^/A (linear 
scale) the system is hardly affected by the relaxation. Between F^/A 
and A (logarithmic scale) the height of current step and conductance 
peak first grow to a maximum value at 6o ~ 0.58F (dashed lines in 
(a) and (b)) for increasing 9, then decrease again and vanish before 
9 — A. Even faster relaxation (reciprocal scale) has no further effect. 
Figures (c) and (d) show cuts through (a) and (b), respectively, for 
five different values of 6. Both the step in (a,c) and peak in (b,d) 
slightly shift towards higher absolute values of eVbias for increasing 
rates between F^/A and A. 



and sub-regimes is responsible for the rather unexpected trans- 
port behavior the quantum dot shows in the presence of spin- 
relaxation. That is to say, the conductance peaks at the onset 
of sequential transport are, as stated above, maximally pro- 
nounced for a small finite relaxation rate. The peaks are situ- 
ated close to the resonances with sequential transitions out of 
the excited state and therefore lie within the intermediate re- 
gion. It turns out that the evolution of the peak height can be 
ascribed to the fact that in the core region of the Coulomb di- 
amond the current I (9) is much more sensitive to changes of 
the relaxation rate 9 than it is in the shell. At small relaxation 
the current is diminished solely in the core region. Hence, the 
height of the current step is increased, while its width remains 
almost constant as compared to zero relaxation. It follows 
that the resulting conductance signatures in the intermediate 
region grow with the relaxation rate as long as the latter stays 
below a level at which the current in the shell region is af- 
fected. 
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At first we describe general features of electron transport 
through the SLQD before we explain in detail how it de- 
pends on the relaxation parameter 9. Fig. 1(b) shows a cal- 
culated charging diagram, i.e., the differential conductance 
G = dl/dVbias against eVbias and e^^^, of the cotunneling 
regime. The parameters are A = Abk-^T, U — 22bk^T, 
r = 4.5 X m-^k^T and 9 = r/2. The Coulomb diamond is 
defined by pronounced red (dark gray) lines of high conduc- 
tance, where one of the reservoir's electro-chemical potentials 
is close to resonance with the energy of a single-charge transi- 
tion involving the ground state. Due to the Coulomb blockade 
the slightly weaker resonance lines of excited state transitions 
are not extended into the diamond. Instead, thin red (dark 
gray) lines appear in that part of the intermediate region, in 
which |eX4)ias| > A. corresponding to the yellow (light gray), 
not dotted area in Fig. 1(a). '"''^ These are the signatures that 
mark the onset of cotunneling-mediated sequential transport 
out of the excited state. The light vertical shades along the 
lines eVbias — ±A arise from the above-mentioned opening 
of inelastic transport channels. Though being actually invalid 
near the border of the Coulomb diamond and in its exterior, 
the approximate solution / was used for Fig. 1(b) even in 
these regions, because its deviation from the systematic solu- 
tion turned out — for the chosen set of parameters — to 
be too small to be visible. In general it may be necessary to 
use the systematic expansion for the border and outer region, 
which can be seamlessly connected to the approximate solu- 
tion in the shell region given that r/3 ^ I''**. 




(a) (b) 

FIG. 3: (Color online). Position of the Fermi levels /xl and /^r of 
left and right reservoir (blue (very dark gray)) relative to the four 
energy differences fioi'P'j'f') ~ ^4>' — ^<f> — e4>_D (light red (light 
gray) horizontal lines) between initial and finial state energy of the 
single-particle transitions (parameters A, U, as in Fig. 1). In (a) 
bold green (light gray) and red (dark grey) arrows represent elastic 
tunneling via virtual state |d) through the dot in ground and excited 
state, respectively, indicated by filled green (light gray) and red (dark 
gray) points (eVbias = 1.5A). Double pointed arrows on the left show 
for both cases how much energy the tunneling electron must gain to 
get into the virtual state, (b) Coherent inelastic processes causing 
transitions |g) |e) (green (light gray) arrows) and \e) ~* |g) (red 
(dark gray) arrows) are illustrated (eVbias ~ 2A). Electrons in the 
red (dark gray) colored (green (light gray) hatched) part of reservoir 
R can tunnel into the dot, if it's in initial state |e) (|g)). 

The dependence of current and conductance on the relax- 
ation rate 9 in the vicinity of the excited state resonance is 
shown in Fig. 2 part (a) and (b), respectively, (the range of 
bias values eVbias is marked by the short horizontal line in Fig. 
1(b)). The parameters A, U, and T are the same as in Fig. 1 



and = {U + A)/2. If the rate 9 is smaller than T^/A, 
thus well below the lowest second-order tunneling rate, the re- 
laxation hardly affects / and G, so that the system behaves as 
for 9 — 0. As the rates grow between F^/A and ~ O.SSF 
they eventually become much larger than every cotunneling 
rate, while still being smaller than sequential rates of energet- 
ically allowed processes. Height and slope of the current step 
as well as the height of the conductance peak increase in this 
range to a maximum value at 9 — 9o, since the relaxation di- 
minishes / in the inelastic part of the core and low bias part 
of the intermediate region but leaves it almost unaltered in 
the sequential part of the shell region. A relaxation with rate 
^^0 < ^ < A is (roughly) as fast as or faster than sequential 
tunneling and while / in this parameter range shows no fur- 
ther relaxation dependence in the core, it decreases in the shell 
region for growing 9. As a result the current step and conduc- 
tance peak decrease as well and vanish before 6' = A ^ F. 
For rates A < 9 < oo transport properties do not depend 
on 9. Obviously there exists an optimal relaxation rate 9a 
for which the signatures of cotunneling-mediated sequential 
transport have maximal height and are considerably more pro- 
nounced than for 9 = 0. Also do the positions of the res- 
onance signatures move to higher absolute values of eVbias, 
when 9 is increased between F^/A and A. As 9o depends on 
the ratio A/U, the coupling F, and temperature /3 in a com- 
plicated way, we can give no simple estimation of its value 
in terms of the system parameters other than: 6*0 = c F with 
0-3 ^ c < 1.2 for all parameter sets yielding reliable results. 
The height of the conductance peaks depends most strongly on 
9 throughout the whole range (F/A)i/2 < 0/y < (A/F)i/2 
with F <C A (see Fig. 2(b)), while its relative variation for all 
values of c amounts to only a few percent. Therefore, we con- 
tent ourself with the statement, that corresponds to a rate 
that is roughly as large as F and very much larger than any co- 
tunneling rate. In order to explain these observations, we ex- 
amine how the tunneling processes that dominate the current 
in the relevant parts of the cotunneling regime are influenced 
by the relaxation. 

In the inelastic part of the core region the current is caused 
solely by elastic and inelastic cotunneling, which also dom- 
inates the occupation of the single-particle states. In par- 
ticular, inelastic tunneling provides an occupation of the ex- 
cited state of order 1 that doesn't depend on F and is re- 
duced by the relaxation, as soon as the magnitude of 9 be- 

(9) 

comes at least comparable to WeV-g. Since the Po.d are much 
smaller than the single-particle occupations in the cotunnel- 
ing regime, the current consists mainly of two contributions 
associated with cotunneling out of state |g) and |e), which 
are proportional to Pg and Pe, respectively. Hence, the co- 
tunneling based on processes with initial state |g) benefits 
from a change dPg{9) « —dPe{9) > caused by relaxation, 
whereas the current with the dot being initially in state |e) 
is decreased by it. So in the core region the dependence of 
the current change dl on dPg{9) is given by dI/dPg{9) « 

—e/h (^^'i^i ~ 1^0''i^e) • sums on the rhs have the 
same sign, which is for both contributions specified by the di- 
rection of current flow and thus by the sign of Hias- For the 
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inelastic part of the core one can establish the relation 



i,(2) 



(13) 



by looking at the energy dependence of elastic and inelas- 
tic processes. In the core, an electron that elastically tunnels 
through the dot with initial (and final) state |x = gj via vir- 
tual intermediate state |d) has to overcome at least the energy 
difference C/ + (5^^gA-e($|5 + |Vbias| /2) > 0, which is byA 
smaller for an initially excited dot than for one in the ground 
state (s. Fig. 3 (a)). The latter is also true for tunneling via 
virtual state |0), which can be seen analogously. As a conse- 
quence the rate for elastic cotunneling is smaller for x = g 
than for x e. 

For the inelastic processes a similar energy argument can be 
applied. Inelastic tunneling out of the ground into the excited 
state cannot set on before e |Vbias| = A, because the energy A, 
needed for the transition to take place, has to be provided by 
the reservoirs. In contrast, inelastic tunneling, causing the op- 
posite transition, is always possible, because in this case the 
transition energy is provided by the dot. Hence, if T4)ias and 
<i>^ specify a point in the core region, for x = e there are 
always more electrons available for inelastic processes com- 
pared to the case x = g (s. Fig. 3 (b)). This results in a higher 
rate for inelastic tunneling through a dot in initial state |e) and 
immediately leads to Eqn. (13). Using this equation and the 
fact that dPg{d) is a positive, monotonic function in 9, for the 
core we can derive 



m 

de 



e dPg 

him 



i,(2) 



< 0. (14) 



Since in this region the single-particle occupations are deter- 
mined mainly by cotunneling, they only depend on the relax- 
ation, if is comparable to the second-order rates. Due to 
the factor dP„/d9 on the rhs of Eqn. (14), this dependence 
also holds for /, whose absolute value decreases for 6 grow- 
ing between F^/A and F and is constant for slower or faster 
relaxation, respectively. 

In the inelastic part of the shell region the maximal 
cotunnehng-provided occupation of the excited state is by 
a factor F smaller than in the inelastic core, since it is re- 
duced by sequential transport out of state |e), as long as 



< W 



L,(l) 



F. When the relaxation becomes faster than 



sequential tunneling, Pg decreases and eventually goes to for 
6* ^ F. The difference in size of sequential and cotunneling 
rates compensates for the reduction of the excited state occu- 
pation, so that the total current is higher in the shell compared 
to the core region. It consists of contributions associated with 
sequential tunneling out of states |x = 0,e,d) and cotunnel- 
ing out of the ground state. 

When 9 is not much larger than F these contributions are 
all of the same order of magnitude, which is, however, not the 
case for their response to increasing relaxation. Obviously, re- 
laxation rates much higher than F completely depopulate the 
excited state and the current caused by tunneling out of |e) 
vanishes. Its relative change in magnitude compared to the 



case of low relaxation is therefore of order 1 . For cotunnel- 
ing out of state |g), on the other hand, the maximum relative 
change is 



Pg(g»F)-Pg(g«F) ^ 0(F) 

pje « r) (1 - o{T)) 



0{T). (15) 



Hence, if the relaxation increases, the gain in the cotunneling 
current associated with Pg cannot compensate for the simul- 
taneous suppression of the sequential current proportional to 
Pg, which results in a reduced total current. Similarly to the 
discussion of the core, however, it can be argued that in the 
shell, where the single-particle occupations are mainly deter- 
mined by sequential processes, the total current can only show 
a considerable relaxation dependence, if 9 is neither much 
smaller nor much larger than F or, equivalently, than the rates 
for sequential tunneling. As we stated above, the current de- 
pendence on 9 smoothly crosses over between the core- and 
shell-like behavior in the intermediate region, so that both the 
current step and the conductance peak grow with 9 between 
F^/A and 9o w cF, while they decrease for 9 > 9o and 
vanish before 9 = A. The fact that the current becomes less 
sensitive to relaxation for higher values of e |X4iias|, showing 
a sharp step-like dependence in the intermediate region, also 
manifests in the slight shift of the position of the excited state 
resonances to higher absolute bias values. 



V. SUMMARY 

In this paper we discussed Coulomb-blocked electron trans- 
port through a SLQD with spin-split level that is coupled 
to two non-magnetic, metallic leads. We used the real-time 
diagrammatic technique to systematically expand occupation 
probabilities and tunneling current up to the second-order in 
the strength F of the tunnel coupling, thereby including se- 
quential and cotunneling into the transport calculations. Two 
properties were considered with respect to which the Coulomb 
blockade regime can be subdivided into parts that differ at 
least in one of them: the kind of tunneling processes con- 
tributing to the current (elastic, inelastic, sequential transport) 
and those determining the single-particle occupations (cotun- 
neling in the core, sequential tunneling in the shell region). 
At or close to the borders between these sub-regimes, signa- 
tures of the dot's excitation spectrum appear in the current 
and differential conductance. With the focus on excited state 
signatures marking the onset of cotunneling-mediated sequen- 
tial transport, we studied how the current is influenced by a 
phenomenologically introduced spin relaxation with rate 9. 
It turned out that for a relaxation rate of about half the tun- 
nel coupling the excited state resonances are maximally pro- 
nounced, being considerably larger than without relaxation, 
while in the limit of infinite 9 the resonances completely van- 
ish. We explained this behavior by a combination of two ef- 
fects: (i) the current decreases monotonically with growing 
relaxation rates and (ii) the excited state occupation is in the 
cotunneling-dominated core and in the shell region only af- 
fected by a relaxation with rates in the range F^/A < 9 < cT 
and cT <9 < A with 0.3 < c < 1.2, respectively. 
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This relaxation dependence of the current may illuminate 
why the resonance signatures measured in Ref. 17 are rel- 
atively sharp compared to the ones that were calculated for 
6 — 0. Furthermore it could provide means to directly influ- 
ence the single-particle occupations in experiments and allows 
to facilitate measurements of excited state resonances by ad- 
justing either the coupling T or the rate 6. 
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